Geometric scaling in inclusive charm production 



V.P. Gongalves x < a and M.V.T. Machado 2 - a - b 
a Institute) de Fisica e Matemdtica, Universidade Federal de Pelotas 

Caixa Postal 354, CEP 96010-090, Pelotas, RS, Brazil 
b High Energy Physics Phenomenology Group, GFPAE, IF-UFRGS 
Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil 

We show that the cross section for inclusive charm production exhibits geometric scaling in a large 
range of photon virtualities. In the HERA kinematic domain the saturation momentum Ql&t{x) stays 
below the hard scale [jq = 4m;?, implying charm production probing mostly the color transparency 
regime and unitarization effects being almost negligible. We derive our results considering two 
saturation models which are able to describe the DESY ep collider HERA data for the proton 
structure function at small values of the Bjorken variable x. A striking feature is the scaling on 
r = Q 2 /QLt(s) above saturation limit, corroborating recent theoretical studies. 
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Introduction. The behavior of ep/pp scattering in the 
limit of high center-of-mass energy ^J~s and fixed mo- 
mentum transfer is one of the outstanding open ques- 
tions in the theory of the strong interactions. Over 
the past few years much theoretical effort has been de- 
voted towards the understanding of the growth of the 
total scattering cross sections with energy. These stud- 
ies are mainly motivated by the violation of the unitar- 
ity (or Froissart) bound by the solutions of the linear 
perturbative DGLAP 1] and BFKL evolution equa- 
tions. Since these evolution equations predict that the 
cross section rises obeying a power law of the energy, vi- 
olating the Froissart bound 3], new dynamical effects 
associated with the unitarity corrections are expected 
to stop its further growth 0, |5j . This expectation can 
be easily understood: while for large momentum trans- 
fer fcj_, the BFKL equation predicts that the mechanism 
.9 ~~ * 99 populates the transverse space with a large num- 
ber of small size gluons per unit of rapidity (the trans- 
verse size of a gluon with momentum k± is proportional 
to l/fcj_), for small fc^ the produced gluons overlap and 
fusion processes, gg — > g, are equally important. Consid- 
ering the latter process, the rise of the gluon distribution 
below a typical scale is reduced, restoring the unitar- 
ity. That typical scale is energy dependent and is called 
saturation scale Q S at- The saturation momentum sets 
the critical transverse size for the unitarization of the 
cross sections. In other words, unitarity is restored by 
including non-linear corrections in the evolution equa- 

tions a a a a i a na m m m q m e 0. 

Such effects are small for k\ > Q 2 at an d verv strong 
for k\ < <5g at , leading to the saturation of the scattering 
amplitude. The successful description of all inclusive and 
diffractive deep inelastic data at the collider HERA by 
saturation models 0, 0, |2(j suggests that these effects 
might become important in the energy regime probed by 
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current colliders. 

One striking feature of the available saturation ap- 
proaches is the prediction of the geometric scaling. 
Namely, the total j*p cross section at large energies is not 
a function of the two independent variables x and Q, but 
is rather a function of the single variable r = Q 2 /Q 2 at - 
As usual, Q 2 is the photon virtuality and x the Bjorken 
variable. A current open question is in what extent the 
geometric scaling is valid above the saturation scale by 
studying the available high energy formulations for the 
linear regime. In Ref. [21| the authors have demon- 
strated that the geometric scaling predicted at low mo- 
menta Q 2 < Q 2 at (x) is preserved by the BFKL evolu- 
tion (at both fixed and running coupling constant) up to 
relatively large virtualities, within the kinematical win- 
dow 1 < lnr ln(Q 2 at /Aq CD ). On the other hand, in 
Ref. [22, the impact of the QCD DGLAP evolution on 
the geometric scaling has been studied. In this case, the 
DGLAP evolution equation is solved imposing as initial 
conditions along the critical line Q 2 — Q 2 at (s) satisfy- 
ing scaling, showing that it is approximately preserved 
at very small x. The residual scaling violation is fac- 
tored out, although the determination of a window for 
the scaling above Q sa t is not provided. As demonstrated 
in Ref. |2lj, the HERA data on the proton structure 
function F2 are consistent with scaling at x < 0.01 and 
Q 2 < 400 GeV 2 . Similar behavior has been observed in 
exclusive [24] and electron- nuclei processes j2^. These 
results, while not entirely compelling, provide a strong 
motivation for further investigations. 

Here we show that the data on inclusive charm pro- 
duction at HERA exhibit geometric scaling above the 
saturation scale. Namely, the charm inclusive cross sec- 
tion depends upon r = Q 2 / Q 2 at (x) alone, where we have 
taken into account that the saturation scale is given by 
the saturation model [lq. i.e. Q 2 at — {xq/x) x . Moreover, 
we extend the symmetric saturation model proposed in 
Ref. |2fj| for the charm case and derive an analytical ex- 
pression for the charm production cross section, which ex- 
plicitly presents geometric scaling. We shown that both 
prescriptions give similar results. 
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Geometric scaling. Lets consider the deep inelastic 
scattering in the dipole frame, in which most of the en- 
ergy is carried by the hadron, while the virtual photon 
7* has just enough energy to dissociate into a quark- 
antiquark pair before the scattering. In this represen- 
tation the probing projectile fluctuates into a quark- 
antiquark pair (a dipole) with transverse separation r± ~ 
1/Q long after the interaction, which then scatters off the 
proton. The interaction 7*p is further factorized in the 
simple formulation [27| . 
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where z is the longitudinal momentum fraction of the 
quark, x ~ Q 2 /W 2 p is equivalent to the Bjorken variable. 
The photon wave functions ^l,t are determined from 
light cone perturbation theory and read as (2*3 
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where the auxiliary variable e 2 = z(l — z) Q 2 + m 2 de- 
pends on the quark mass, rrif. The Kqi are the McDon- 
ald functions and the summation is performed over the 
quark flavors. 

The dipole hadron cross section a d ip contains all in- 
formation about the target and the strong interaction 
physics. In the Color Glass Condensate (CGC) formal- 
ism 0, [n| , Udip can be computed in the eikonal approx- 
imation, resulting 



Pdip(x 



.) = 2 J d 2 b ± [l-S(z,ru,M] . 



(2) 



where S is the ^-matrix element which encodes all the in- 
formation about the hadronic scattering, and thus about 
the non-linear and quantum effects in the hadron wave 
function. The function S can be obtained by solv- 
ing an appropriate evolution equation in the rapidity 
y = ln(l/a;). The main properties of S are: (a) for the 
interaction of a small dipole (r± <C 1/Qsat), S(r±) « 1, 
which characterizes that this system is weakly interact- 
ing; (b) for a large dipole (r± 3> l/Q S at), the system is 
strongly absorbed which implies S(r±) <C 1. This prop- 
erty is associate to the large density of saturated gluons 
in the hadron wave function. In our analysis we will 
initially consider the phenomenological saturation model 
proposed in Ref. |18| which encodes the main properties 
of the saturation approaches. In this model 
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with o~dipl o~o the scattering amplitude, averaged over all 
impact parameters b±, and Q 2 at ~ A 2 e xin ( x o/x) _ The pa- 
rameters of the model were constrained from the HERA 



small x data, coming out typical values of order 1-2 
GeV 2 for the momentum scale. We have that when 
Q 2 at (a;) r\ <C 1, the model reduces to color transparency, 
whereas as one approaches the region Q 2 at (a;)r 2 L w 1, 
the exponential takes care of resumming many gluon ex- 
changes, in a Glauber-inspired way. Intuitively, this is 
what happens when the proton starts to look dark. It 
is important to emphasize that one moves towards the 
unitarity bound for large r± in this saturation model 
much faster than the predicted by the CGC approach. 
However, as our goal is the charm production, which is 
dominated by the color transparency regime, the differ- 
ence among the approaches can be disregarded in what 
follows. The saturation model depends upon the vari- 
ables x and r± only through the dimensionless quantity 
Qsat T± . Consequently, the saturation model predicts the 
geometric scaling of the total cross section. 

We are interested in the charm production in deep in- 
elastic scattering. From the experimental point of view, 
the HERA experiments have published data for the con- 
tribution of charmed meson production to the structure 
function F%. Their analysis were based on D° and D* 
meson tagging. This allows one to single out the charm 
contribution i 7 ^ to the total structure function and thus 
to investigate if the property of geometric scaling is also 
present in this observable. Before presenting our results, 
lets perform a qualitative analysis of the inclusive charm 
production using the saturation model |l8j in order to 
shed light on the dipole configurations dominating the 
process in the relevant kinematical limits and show how 
the geometric scaling comes out. A characteristic fea- 
ture in heavy quark production within the color dipole 
approach is that the pro cess is dominated by small size 
dipole configurations [28| . The overlap function weight- 
ing the dipole cross section is peaked at r± ~ l/m c ~ 0.1 
fm even for sufficiently low Q 2 values. As a conse- 
quence, charm production is dominated by color trans- 
parency and saturation effects are not important there, 
i.e. cr dip ~ S (x,r ± ) ee (T Q 2 at (a;)r 2 L /4, which is the 
leading order term in the expansion of S for small r±. 
We consider for simplicity the transverse contribution, 
whose corresponding wave function is given by Eq. ffl. 
and taking its asymptotic behavior at e r± <C 1 such that 
K±(sr±) ~ l/sr±. This imposes additional constraint to 
the z integration, namely we have the integrand should 
be multiplied by the Heaviside function 0(1— e 2 rj_), then 
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where fi 2 — Am 2 . The first term corresponds to the sym- 
metric dipole configurations, i. e. < z >« 1/2, whereas 
the second term comes from aligned jet configurations, 
with the charm mass introducing a cut-off on the maxi- 
mal size of the charmed dipole. 

For the HERA kinematical region, we have Q 2 at w 1 
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FIG. 1: Experimental data on inclusive charm production 
plotted versus the scaling variable r = Q 2 /<2sat- 



FIG. 2: The quantity ^/t a to t plotted versus the scaling vari- 
able r. 



GeV 2 , which implies that the relation Q 2 at < Q 2 + y? c 
is ever satisfied. Consequently, we can define two kine- 
matical regimes depending of the relation between Q 2 
and fi 2 . For Q 2 ^> /i 2 we have scaling with logarith- 
mic enhancement coming from aligned jet configurations, 
whereas for Q 2 <C /x 2 only symmetric dipole configura- 
tions contribute. Therefore, we obtain that 
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where the first term provides the behavior 1/r at large 
r whereas the second term leads to a smooth transi- 
tion down to the asymptotic (r-independent) behavior 
at small t. 

An analytical expression for the r dependence of the 
inclusive charm production can also be obtained in a less 
model dependent way. For this purpose we will make 
use of the symmetric saturation model |26j |. where the 
energy evolution of the proton leads to the parton multi- 
plication and the transverse momentum scale Q S at (x) ap- 
pears. The main assumption is that the evolved proton 
can be described by a collection of independent dipoles 
at the time of the interaction whose sizes are distributed 
around 1 / Q sat . The rate of growth of the parton densities 
is assumed to be Q 2 at (x) /A 2 and the symmetry between 
low and high virtualities in jp interactions comes from 
the symmetry in the dipole-dipole cross section. The ap- 
proach provides a quite intuitive and simplified expres- 
sion for the inclusive production. Following such an ap- 
proach, we obtain that in the HERA kinematical regime, 
i.e. Q 2 at < (jl 2 , the inclusive charm production is given 
by, 
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where r c = fi 2 /Q 2 at (x) and the parameters are taken 
from the data fit in Ref . 26] . Here, we make the simpli- 
fied assumption that the coupling with the dipole is flavor 



blind, in such way that N c£ = (2/5)N, with N being the 
global normalization describing F^ data. The factor 2/5 
corresponds to the charge fraction e 2 /(y^e 2 d s + e 2 ). 
Once the parameters are fitted to proton structure func- 
tion data, our prediction for the r dependence in the 
inclusive charm production is parameter free. Moreover, 
the expression above gives an analytical result for the 
whole r domain. 

Results. In Fig. ^we show the experimental data [2jJ 
on the total cross section for the inclusive charm produc- 
tion plotted versus scaling variable r = Q 2 /Q 2 at , with 
Q S at from the saturation model. We include all available 
data covering the range 1.5 < Q 2 < 130 GeV 2 [2<|. It 
should be stressed that some bins with x > 10~ 2 have 
also been included and old measurements contain some- 
what large uncertainties. Moreover, charm data are in- 
fluenced by a significant systematic uncertainty coming 
from the theoretical models used in the extraction of F^ c . 
The combination of these factors produces a larger data 
dispersion than in the inclusive structure function. We 
see the data exhibit geometric scaling for the whole Q 2 
range, verifying a transition in the behavior on r of the 
cross section from a smooth dependence at small r and 
an approximated 1/t behavior at large r. The transi- 
tion point is placed at ^i 2 = 4m 2 , which takes values of 
order 10 GeV 2 for a charm mass m c — 1.5 GeV. This 
turns out in r ~ 10 since at HERA Q 2 at ~ 1 GeV 2 . The 
asymptotic 1/r dependence reflects the fact the charm 
production cross section scales as Q 2 at /Q 2 modulo a loga- 
rithmic correction ~ ln(<5 2 //z 2 ), with energy dependence 
driven by the saturation scale. The mild dependence at 
t < corresponds to the fact the cross section scales 
as Q 2 at //i 2 towards the photoproduction limit, but with 
the same energy behavior given by the saturation scale. 
As also plotted in Fig. |2 we also found a symmetry 
between the regions of large and small r for the func- 
tion y/r <J% t with respect the transformation r «-* 1/r 
in the whole region of r. The features present in the 
inclusive charm production data can be well reproduced 
in the phenomenological saturation model as shown in 



4 



Eq. J5J, corresponding to the solid curve in Figs. ^ an d 
13 The symmetric saturation model also provides similar 
results, as shown in the dot-dashed lines. Disregarding 
the Glauber-like resummation in this model, the expres- 
sion gets simplified to af ot cx t ^_ t [1 + log(r + r c )], and 
the symmetric pattern is easily verified. Moreover, it 
is important to emphasize that a reasonable description 
of the charm photoproduction experimental data, which 
corresponds to r = 0, is also obtained using the symmet- 
ric saturation model. Our results for the nuclear heavy 
quark photoproduction indicate that similar behavior is 
expected in the nuclear case [30|. 

Summary. Summarizing our results, one shows that 
the inclusive charm production exhibits geometric scal- 
ing in a large region of photon virtualities. In the HERA 
kinematic domain the saturation momentum Q 2 &t (x) 
stays below the hard scale [± 2 — 4m;?, implying that 
charm production probes mostly the color transparency 
regime and saturation corrections are not very important. 



In the color dipole picture, the transition at Q 2 ~ fx 2 is 
related to the presence of aligned jet configurations at 
Q 2 3> fJ. 2 and complete dominance of symmetric con- 
figurations at Q 2 <C n 2 - An analytical result on the r 
behavior was obtained within the symmetric saturation 
model, relying on quite simple assumptions about the 
dipole-proton interaction. An outstanding feature is the 
scaling on r above saturation scale, supported by recent 
theoretical formulations. 
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